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Abstract 

Haplodiploid sex determination allows unmated females to produce sons. Con- 
sequently, a scarcity of males may lead to a significant proportion of females 
remaining unmated, which may in turn give rise to a surfeit of males in the fol- 
lowing generation. Stable oscillation of the sex ratio has been predicted by clas- 
sic models, and it remains a puzzle as to why this is not observed in natural 
populations. Here, I investigate the dynamics of sex allocation over ecological 
and evolutionary timescales to assess the potential for sustained oscillation. I 
find that, whilst stable oscillation of the sex ratio is possible, the scope for such 
dynamical behavior is reduced if sex allocation strategies are evolutionary labile, 
especially if mated females may facultatively adjust their sex allocation accord- 
ing to the present availability of mating partners. My model, taken together 
with empirical estimates of female unmatedness in haplodiploid taxa, suggests 
that sustained oscillation of the sex ratio is implausible in natural populations. 
However, this phenomenon may be relevant to artificially introduced biological 
control agents. 
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Introduction 

Around 20% of all animal species employ haplodiploid sex 
determination (Crozier and Pamilo 1996). In such species, 
whilst females are produced sexually, deriving from eggs 
that have been fertilized by sperm, males are produced asex- 
ually, deriving from unfertilized eggs, in a process termed 
"arrhenotoky". Thus, haplodiploidy allows unmated 
females to produce offspring, of whom all are male. 

Nearly a century ago, the insect ecologist C. B. Williams 
(1917) highlighted the potential for arrhenotoky to drive 
oscillation of the sex ratio. Specifically, a scarcity of males 
in one generation may lead to a significant proportion of 
females being unmated, which in turn leads to a surfeit of 
males in the next generation. Fifty years later, W. D. Ham- 
ilton (1967) provided an analytical treatment, reporting 
that this would lead to sustained oscillation of the sex ratio, 
provided k < (1— e)/e where k is the number of females 
than can be mated by each male and e is the sex ratio 
(proportion male) among the offspring of mated females. 
Such sustained oscillation obtains even when males are 



highly fecund (large k), provided that the sex allocation of 
mated females is sufficiently female biased (low e) . 

It remains a puzzle as to why such sustained oscillations 
are not, in fact, observed in the natural world. One possibil- 
ity is that Williams' and Hamilton's predictions are artifacts 
of artificial model assumptions. For example, both authors 
treated the sex allocation of a mated female as a fixed 
parameter and did not consider its evolutionary dynamics. 
However, sex allocation does evolve by natural selection, 
and the presence of unmated females who are constrained 
to produce only sons favors mated females to allocate more 
resources to daughters (Godfray and Grafen 1988; Godfray 
1990; West 2009). Consequently, the sex allocation of 
mated females is expected to co-evolve with the extent of 
female unmatedness, making some regions of Hamilton's 
parameter space more evolutionarily plausible than others. 
Thus, the true scope for female unmatedness to drive stable 
oscillation of the sex ratio under haplodiploidy is unclear. 

Here, I extend Williams' and Hamilton's analyses by 
investigating evolutionary change in sex allocation and its 
impact on the ecological dynamics of the sex ratio and 
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female unmatedness. First, I consider that mated females 
employ a fixed sex allocation strategy that is evolution- 
arily optimized according to the intergenerational average 
availability of mating partners. This is consistent with 
Hamilton's ecological model and clarifies the parameter 
range over which stable oscillation obtains. Second, I con- 
sider that mated females are able to facultatively adjust 
their sex allocation according to current availability of 
mating partners. This requires an extension of Hamilton's 
ecological model and further clarifies the parameter range 
over which stable oscillation obtains. 

Models and Results 
Obligate sex allocation 

Following Hamilton (1967), I consider an infinite, haplo- 
diploid population. Mating occurs as follows: females and 
males mate at random, with females leaving the mating 
pool upon their first mating and males leaving the mating 
pool upon their A:th mating, such that mating ends when 
one or both of the sexes is no longer present in the mating 
pool. Mated females produce offspring with sex ratio e, and 
unmated females produce offspring with sex ratio 1, with 
all females producing the same number of offspring. 

Thus, if the sex ratio in any generation is z, the pro- 
portion of females being mated in that generation is 
m = (zfc)/(l-z) if z < 1/(1 + it) and m = 1 if z > 1/ 
(1 + k), and the sex ratio in the following generation is 
zf = m e + l — m, or: 



l-^(l-e) z< 
if 



z > 



i 

1 t-k 
1 

1 I k 



(1) 



As Hamilton reported, this leads to stable oscillation in 
sex ratio between z = e and z = 1— ek, either side of an 
unstable equilibrium at z* = !4(2 + k(l— e)— (fc(l— e) 
(4 + k(l-e)))' A ), if k < (l-e)/e (see Appendix Al for 
derivation; Fig. 1). However, the sex allocation of mated 
females is not arbitrary. Rather, this trait is expected to 
be under strong selection to counter the male bias intro- 
duced by the reproduction of unmated females (Godfray 
and Grafen 1988; Godfray 1990; West 2009). If mated 
females adopt a sex allocation e that is not facultatively 
adjusted according to mate availability, then natural selec- 
tion leads this to converge upon: 



' j + k-y/(9-k)(l-k) 



8k 



k< 1 



k> 1 



(2) 



Stable equilibrium 



Stable oscillation 



(see Appendix Al for derivation; Fig. 2A). This result 
extends that of Godfray (1990), who assumed a fixed rate 
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Male fecundity (k) 

Figure 1. Hamilton's model. Stable oscillation of the sex ratio 
between z = e and z = 1 -ek obtains when k < (1 -e)/e, and a stable 
equilibrium at z* = '/ 2 (2 + k(\ -e)-(k(1 -e)(4 + fc(1 -e))) v ') obtains 
when k > (1-e)/e. Insets illustrate the scenarios indicated by arrows 
(k = 0.10 & e = 0.50 and k = 1.50 & e = 0.75). 

of female unmatedness, to a scenario in which unmated- 
ness fluctuates over generations. Whilst Godfray (1990) 
showed that mated females are favored to adopt a female- 
biased sex allocation that exactly negates the male bias 
introduced by female unmatedness, here, I find that 
mated females are forced to overcompensate when male 
fecundity is limiting and undercompensate when male 
fecundity is not limiting. 

This development is compatible with Hamilton's model 
of the ecological dynamics of the sex ratio, but reduces its 
parameterization. In particular, whereas Hamilton's model 
is governed by two parameters - male fecundity (k) and 
mated-female sex allocation (e) - 1 have expressed the latter 
in terms of the former, so that the model is solely governed 
by the male-fecundity (k) parameter. Returning to Hamil- 
ton's original condition k < (1— e)le for stable oscillation, 
and substituting e for the solution given in equation (2), 
the condition for stable oscillation becomes k < 1. In other 
words, stable oscillation obtains only if each male can, on 
average, mate with fewer than one female. In this case, the 
sex ratio stably oscillates between z = (3 + k— ((9— k) 
(l-k)) y ')/{8k) and z = (5-k+((9-fc)(l-fc))' /j )/8; that is, 
the sex ratio is alternately female biased and male biased in 
successive generations, and on average, it is male biased 
(Fig. 2B). Otherwise, if k > 1, all females are guaranteed of 
being mated when the sex ratio is z = 'A, and the adoption 
of the sex allocation strategy e = l A ensures that the sex 
ratio z* = V2 remains stable (Fig. 2B; cf Williams 1917). 

Facultative sex allocation 

I now assume that mated females produce offspring with 
sex ratio e(m), which is facultatively adjusted according 
to the proportion m of females that are mated in their 
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Figure 2. Obligate sex allocation. (A) Natural selection favors mated females to exhibit female-biased sex allocation e = (3 + k-((9-k){ 1-k)) Y ')/ 
(8k) when k < 1 and unbiased sex allocation e = Vi when k > 1. (B) This leads to stable oscillation of the sex ratio between z = (3 + k-((9-k) 
(1~k)) y ')/(8k) and z = (5-k + {(9~k)( 1-k))' /2 )/8 when k < 1, and a stable equilibrium at z* = i4 when k > 1. Insets illustrate the scenarios 
indicated by arrows = 0.50 and <r = 1.50). 



generation. Natural selection favors the sex allocation 
strategy: 



0 



m < 



e(m) 



(3) 



2m 1 

2m 



(see Appendix Al for details; Fig. 3A). This is Godfray's 
(1990) result, and here, I have shown that it extends to a 
more complicated scenario than the one he considered. 
Specifically, whilst Godfray (1990) assumed that the pro- 
portion of mated females remains fixed over successive 
generations, I have considered that it may fluctuate 
between generations. Crucially, I have considered that 
females respond facultatively to this fluctuation, rather 
than adopting a fixed sex allocation that optimizes with 
respect to average mate availability (as was done in the 
previous section). 

This development is not compatible with Hamilton's 
model of the ecological dynamics of the sex ratio, which 
assumes that mated females adopt the same sex allocation 
in all generations, irrespective of the current availability of 
males. Consequently, I extend Hamilton's model to con- 
sider intergeneration variation in the sex allocation of 
mated females. As in Hamilton's model, if the sex ratio in 
any generation is z, then the proportion of females being 
mated in that generation is m = (zk)l(l-z) if z < 1/ 
(1 + k) and m = 1 if z > 1/(1 + k), and the sex ratio in 
the following generation is z' = m e{m) + 1—m. Making 
the substitution given in equation (3) yields: 



1 



zk 
1-z 



Z< 



z > 



1 

1+2/c 



1 

1+2*: 



(4) 



This leads to stable oscillation in sex ratio when k < 'A. 
In this case, the sex ratio stably oscillates between z = 'A 
and z = l—k; that is, the sex ratio is alternately unbiased 
and male biased in successive generations, and on aver- 
age, it is male biased (Fig. 3B). Otherwise, if k > 'A, 
mated females are guaranteed to be able to set the popu- 
lation sex ratio in the next generation to z = 'A given that 
it is currently at z = 'A, and hence, the sex ratio z* = l A 
remains stable (Fig. 3B). 

Discussion 

Building upon the numerical model of Williams (1917), 
Hamilton (1967) showed that female unmatedness could 
drive sustained ecological oscillation of sex ratio in haplo- 
diploid populations. Here, I have assessed the impact that 
evolutionary optimization of the sex allocation behavior of 
mated females has upon such dynamics. Further to Hamil- 
ton's suggestion that sustained oscillation obtains even for 
arbitrarily large male fecundity, so long as mated-female 
sex allocation is sufficiently female biased, I have found 
that evolutionary optimization of mated-female sex alloca- 
tion behavior places strict upper limits upon male fecun- 
dity, above which stable oscillation does not obtain. 

First, I considered that mated females optimize their 
sex allocation behavior according to the average level of 
female unmatedness experienced over evolutionary time- 
scales (obligate sex allocation). I found that stable oscil- 
lation in the sex ratio obtains only when each male can, 
on average, mate with fewer than one female. Otherwise, 
an even sex ratio guarantees that all females are mated, 
and this ensures that all females - adopting an even sex 
allocation strategy - can recover an even sex ratio in the 



) 201 4 The Author. Ecology and Evolution published by John Wiley & Sons Ltd. 



1625 



Sex Ratio Dynamics under Haplodiploidy 



A. Gardner 




Proportion of females mated (m) Male fecundity (k) 



Figure 3. Facultative sex allocation. (A) Natural selection favors mated females to exhibit female-biased sex allocation e(m) = 0 when m < '/•> and 
e(m) = (2 m-1)/(2 m) when m > V 2 , where m is the proportion of females that are mated in this generation. (B) This leads to stable oscillation of 
the sex ratio between z = 'A and z = 1 -k when k < V%, and a stable equilibrium at z* = !/ 2 when k > V%. Insets illustrate the scenarios indicated by 
arrows (k = 0.25 and k = 1.50). 



subsequent generation. In the event that male fecundity 
is sufficiently low for stable oscillation in the sex ratio 
to obtain, mated females are favored to adopt a female- 
biased sex allocation strategy that undercompensates for 
the male bias introduced by unmated females in those 
generations in which there is a scarcity of males, and 
overcompensates in those generations in which there is a 
surfeit of males. The result is a sex ratio that oscillates 
between male bias and female bias in alternate genera- 
tions, with an average male bias. 

Second, I considered that mated females optimize 
their sex allocation behavior according to the present 
level of female unmatedness experienced in their genera- 
tion (facultative sex allocation). I found that stable oscil- 
lation in the sex ratio obtains only when each male can, 
on average, mate with fewer than 0.5 females. Otherwise, 
an even sex ratio guarantees that at least half of all 
females are mated, and this ensures that all females - 
adopting a suitably female-biased sex allocation strategy 
- can recover an even sex ratio in the subsequent gener- 
ation. In the event that male fecundity is sufficiently low 
for stable oscillation in the sex ratio to obtain, mated 
females are favored to adopt a female-biased sex alloca- 
tion strategy that exactly compensates for any male bias 
owing to female unmatedness. 

My obligate sex allocation analysis generalizes that of 
Godfray (1990), who considered optimization of mated- 
female sex allocation in a constant mating environment, 
without sex- ratio oscillation. Godfray (1990) found that 
mated females are favored to adopt a sex allocation of 
e = (2 m— 1)/(2 m), where m is the frequency of mated 
females. In the context of fluctuations in rates of unmat- 
edness, this becomes e « ((2m — l)/(2m )) + (a,„ 2 l(2m 3 )) 



where m is the average and er m 2 the variance in the fre- 
quency of mated females (see Appendix Al for details). 
My facultative sex allocation analysis yields a result 
identical to that given by Godfray (1990), but in a novel 
context. Specifically, if mated females employ obligate sex 
allocation in a constant mating environment (Godfray 
1990), or facultative sex allocation in a variable mating 
environment (this paper), they are favored to adopt the 
strategy e = (2 m— 1)/(2 m). 

These developments build extra realism into Hamil- 
ton's (1967) model of sex- ratio dynamics. Usually, greater 
realism is achieved by adding parameters to models, mak- 
ing them more complex. However, greater realism has 
been achieved here by removing a parameter, by consider- 
ing that sex allocation is an evolving variable. Other ways 
in which the model could be made more realistic include 
considering that the number of matings per male may be 
dependent upon the sex ratio, rather than taking a con- 
stant value, and allowing females imperfect information 
as to their mating environment. I leave these as avenues 
for future exploration. 

Evolutionary optimization of sex allocation behavior 
greatly restricts the conditions under which stable oscil- 
lation may occur in natural populations. If mated 
females are obliged to employ a fixed sex allocation 
behavior, then even a small proportion of females 
remaining unmated may be sufficient to drive sustained 
oscillation in the sex ratio. But if mated females can fac- 
ultatively adjust their sex allocation according to current 
male availability, then sustained oscillation requires that 
a majority of females remain unmated in every genera- 
tion. Empirical estimates of the incidence of unmated 
(or otherwise "constrained") females in haplodiploid 
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taxa range from 0 to 30% and are usually <5% (West 
2009), so there appears to be very little scope for sus- 
tained sex-ratio oscillation in natural populations. How- 
ever, unmated females could plausibly outnumber mated 
females in artificially introduced biological control spe- 
cies (Rhainds 2010), especially at range frontiers and if 
dispersal is female biased (Heimpel and Asplen 2011). 
Haplodiploid taxa - especially parasitoid wasps - have 
been particularly favored for controlling agricultural 
pests (Waage and Hassell 1982). More generally, even 
transient oscillation may be long lasting over economi- 
cally relevant timescales, making this an important factor 
in the design of biological control programs. 
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Appendix A1: Dynamics of sex ratio 
and female unmatedness 

Obligate sex allocation 

Ecological dynamics 

Inspection of equation ( 1 ) suggests a candidate equilibrium 
at z* = e, and this requires e > 1/(1 + k). If e + 8 > II 
(1 + k), where 8 is a vanishingly small quantity, then z 1 = e 
when z = e + 5, so the equilibrium z* = e is stable if 
e > 1/(1 + k) or, equivalently, k > (1— e)le. A second can- 
didate equilibrium is found by setting z! = 1— (zfc/(l— z)) 
(1— e) = z and solving for z = z* to yield z* = Vz{2 + k 
(l-e)-(fc(l-e)(4 + k{\-e)))' A ). This assumes z* < 1/ 
(1 + k), which is equivalent to k < (1— e)le. Setting z = z* 
+ <5 yields z' = z*— (l/z*)<5 + 0(<5 2 ) which, as the magni- 
tude of coefficient of 8 is never less than unity, means that 
the equilibrium is unstable. Finally, if the sex ratio in any 
generation is z = e, then the proportion of females that are 
mated is ek/(l—e), and the sex ratio in the subsequent gen- 
eration is zf = l—ek, which is sufficient for all females in 
that generation to be mated because k < (l—e)/e. Hence, 
the sex ratio in the generation after that is z" = e. Thus, the 
sex ratio oscillates between e and l—ek. 

Evolutionary dynamics 

Sex allocation can be modeled by assigning every female 
an equal number of sons and daughters and then elimi- 
nating half of the offspring of each female according to 
her mating status and her sex allocation strategy (Taylor 
and Frank 1996). Thus, the probability that a focal female 
survives this culling is m(l— x), and the probability that a 
focal male survives this culling is mx + 1— m, where m is 
the probability that their mother is mated and x is the 
proportion of their mother's reproductive resources that 
is invested into sons in the event that she is mated. The 
population average female survival is m(l— e), and the 
population average male survival is me + 1— m. Hence, 
the relative fitness of a focal female is W t = (m(l—x))/(m 
(1— e)) = (1— x)l(l— e), and the relative fitness of a focal 
male is W m = {mx + l—m)/(me + l—m). Natural selec- 
tion maximizes W = £ m n„,(c{ W f + c m W m ), where n m is 
the proportion of generations in which the frequency of 
mated females is m, and c f = 2/3 and c m = 1/3 are the 
class reproductive values of females and males, respec- 
tively (Fisher 1930; Taylor 1996; Taylor and Frank 1996). 
That is, if a focal locus controls the sex allocation e of 
mated females, then, drawing a gene from this locus at 
random from the population and denoting its genie value 
by g, the condition for natural selection to favor an 
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increase in e is dW/dg > 0, that is, £m n m (c( dWf/dgf + c m 
dWJdg m ) > 0. Here, dW f /% = {dWf/dx) x (dxl 
dG) x (dG/dgf), where G is the genetic "breeding" value 
of the individual's mother for the sex allocation pheno- 
type, dx/dG = 1 is the genotype-phenotype map and dG/ 
dgf = % is the consanguinity of mother and daughter 
(Taylor and Frank 1996). Similarly, dW m /dg m = (9W m / 
ck) x (dx/dG) x (dG/dg m ), where dG/dg m = Vz is the 
consanguinity of mother and son. This yields the condi- 
tion Y,m n m (m/(m e + l—m)) - (1/(1— e)) > 0. In gen- 
eral, this leads to a convergence stable value e w ((2 m 
— l)/(2m )) + (cr m 2 /(2m 3 )), where m is the average and 
o m 2 the variance in the proportion of females being 
mated (the approximation neglects higher orders of 
variance). However, of immediate interest is the scenario 
where the sex ratio oscillates between e and 1— ek, so that 
m takes values mi = ek/(l—e) and m 2 = 1 with equal 
probability. Thus, the condition for increase in e is (m,/ 
(wii e + l-m-J) - (l/(l-e)) + (m 2 /(m 2 e + l-m 2 )) - (1/ 
(1— e)) > 0, which yields an exact convergence stable 
value for e, given by equation (2). 

Facultative sex allocation 
Evolutionary dynamics 

The expected fitness of a juvenile female whose mother's 
generation involved a proportion m* of females being 
mated is W{ = m* (1— x), and the expected fitness of a 
juvenile male whose mother's generation experienced a 
proportion m* of females being mated is w m = m* 
x + l—m*, where x is the proportion of their mother's 
reproductive resources that is invested into sons in the 
event that she was mated. Noting that the population aver- 
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age of x is x = e(m*), the average fitness for females and 
males in such generations is Wf = m*(l— e(m*)) and 
w m = m* e(m*) + l—m*, respectively. As before, the con- 
dition for natural selection to favor an increase in e(m*) is 
X„, n rn (c{ dW { /dg { + c m dWJdg m ) > 0, and as the gene is 
only active when m = m*, this reduces to n m * (cf dWf/ 
dgf + c m dW m /dg m ) > 0 where derivatives are evaluated at 
m = m* and x = x. This yields the condition -(1/(1— e 
(m*))) + (m*/(m* e(m*)+l—m*)) > 0, such that the con- 
vergence stable sex allocation strategy for mated females is 
given by equation (3). 

Ecological dynamics 

Inspection of equation (4) suggests the candidate equilib- 
rium z* = 14, and this requires Vz > 1/(1 + 2k), that is, 
k > I/2. If Vz + S > 1/(1 + k), where 8 is a vanishingly 
small quantity, then z! = 'A when z = Vz + d, so the equi- 
librium z* = 'A is stable if k > Vz. A second candidate 
equilibrium is found by setting z! = I— (zkl(l-z)) = z, 
and solving for z = z* to yield z* = 1/2(2 + k— ((4 + k) 
k)' A ). This assumes z* < 1/(1 + 2k), which is equivalent 
to k < 'A. Setting z = z* + 8 yields z 1 = z*— (4/ 
(fc' /2 -(4 + ic)' /2 ) 2 )<5 + 0(<5 2 ) which, as the magnitude of 
the coefficient of 8 is never less than unity, means that 
the equilibrium is unstable. Finally, if the sex ratio in any 
generation is z = Vz, the proportion of females that are 
mated is m = k, and hence, the sex ratio in the subse- 
quent generation is z! = m e(m) + 1— m = I— k. This 
leads to a proportion m = {(l—k)lk)k = l—k of females 
being mated, and consequently, the sex ratio in the subse- 
quent generation is z" = m e(m) + l—m = Vz. Thus, the 
sex ratio oscillates between Vz and l—k. 
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